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Abstract 

In light of recent remarkable progress in Higgs search at the LHC, we study the rare decay 
process h — > and show its correlation with the decay h — > 77 in low energy SUSY models such 
as CMSSM, MSSM, NMSSM and nMSSM. Under various experimental constraints, we scan the 
parameter space of each model, and present in the allowed parameter space the SUSY predictions 
on the Z'j and 77 signal rates in the Higgs production at the LHC and future e + e~ linear colliders. 
We have following observations: (i) Compared with the SM prediction, the Zj and 77 signal rates 
in the CMSSM are both slightly suppressed; (ii) In the MSSM, both the Z7 and 77 rates can be 
either enhanced or suppressed, and in optimal case, the enhancement factors can reach 1.2 and 2 
respectively; (iii) In the NMSSM, the Z7 and 77 signal rates normalized by their SM predictions 
are strongly correlated, and vary from 0.2 to 2; (iv) In the nMSSM, the Z7 and 77 rates are 
greatly reduced. Since the correlation behavior between the Z7 signal and the 77 signal is so 
model-dependent, it may be used to distinguish the models in future experiments. 

PACS numbers: 14.80.Da,14.80.Ly,12.60.Jv 
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I. INTRODUCTION 



Based on the measurements of 77 and ZZ* channels the ATLAS and CMS collaborations 
have independently provided compelling evidence for a bosonic resonance around 125-126 
GeV jl|. This is a great triumph for particle physics, but it also leads to a host of new 
questions about the nature of the boson. So far there exist large uncertainties in determining 
the rates of the two channels, and meanwhile, observation of the boson through other signals, 
such as bb and r + r~ channels, is still far away from becoming significant j^j] . So although the 
preliminary data of the LHC indicate that the boson closely resembles the Higgs boson in 
the Standard Model (SM), the deficiencies of the SM itself such as gauge hierarchy problem 
suggest new physics explanation of the boson. Obviously, in order to decide the right 
underlying theory, the LHC should exhaust its potential to measure the decay channels of 
the boson as accurately as possible in its high luminosity phase. 

Among the decay modes of the Higgs-like boson h, the diphoton channel plays a very 
important role in determining its mass, spin and parity ^|. At the same time, since the 
diphoton channel is mediated by loops of charged particles, it also acts as a sensitive probe 
to new physics. In fact, this feature has been widely utilized to explain the excess of the 
diphoton signal over its SM Higgs boson hypothesis reported by the ATLAS and CMS 
collaborations (4 -14]. In this note, we concentrate on another decay mode h — > Z7. In the 
SM, the branching ratio of this decay is about two thirds of that for the diphoton decay, and 
just like the diphoton signal, it can provide a clean final-state topology in determining the 
properties of the boson, such as its mass, spin and parity jl5|. Moreover, since new charged 
particles affecting the diphoton decay can also contribute to the Z7 decay, the two decay 
modes should be correlated, and therefore studying them in a joint way can reveal more 
details about the underlying physics. Albeit the advantages, in contrast to the diphoton 
decay which has been intensively studied, the Z7 decay was paid little attention in the past, 
or example, so far only a few works have been devoted to this decay in new physics models 
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171 ] . and until very recently have the CMS and ATLAS collaboration set an upper limit 
on the ratio Czy/v^y < 10[18j. Here we remind that, although the Z7 signal suffers from 
large irreducible background at the LHC 18|], at the next generation linear colliders with 



given collision energy, the Higgs event from the process e + e — > Zh — > ZZj can be easily 
reconstructed, which is very helpful in suppressing the background for such signal. So there 
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is good prospect to measure accurately this decay in future. 

In the following, we focus on the Z7 decay channel of the SM-like Higgs boson h in low 
energy supersymmetric models such as the Constrained Mimimal Supersymmetric Stan- 
dard Model (CMSSM), the Minimal Supersymmetric Standard Model (MSSM), the Next- 
to-Minimal Supersymmetric Standard Model (NMSSM) and the Nearly Minimal Supersym- 
metric Standard Model (nMSSM). We investigate the Z7 signal of the Higgs production at 
the LHC and future e + e~ linear colliders, and especially, we show its correlations with the 
77 signal. As we will show below, the Z7 signal rate may be either enhanced or suppressed 
in SUSY, and its correlation behavior is so model-dependent that it may be utilized to 
distinguish the models in high luminosity phase of the LHC. 

This work is organized as follows. In Section II we introduce the basic features of the 
SUSY models and present some formulae relevant to our calculation. In Section III we first 
discuss the effects of new charged SUSY particles on the partial decay widths of h — > Zj, 
then we study in a comparative way the Z7 and 77 signal rates of the Higgs production at 
different colliders. Finally, we draw the conclusions in Section IV. Various couplings used in 
the calculation are given in the Appendix. 

II. THE MODELS AND ANALYTIC FORMULAE 

In a low energy supersymmetric gauge theory, the explicit form of its Lagrangian is 
determined by the gauge symmetry, superpotential and also soft breaking terms. As for the 
four models considered in this work, their differences mainly come from the superpotential, 
which is the source for the Yukawa interactions of fermions and self interactions of scalars. 

MSSM and CMSSM : The MSSM and CMSSM contain two Higgs doublets H u , H d 
and so they predict five physical Higgs bosons, of which two are CP-even, one is CP-odd 
and two are charged. Their superpotential is given by 

W MSSM = W F + ^H U -H d , (1) 

where Wf denotes the Yukawa interaction, and it takes following from 

W F = 1Y U Q ■ H u - ~dY d Q ■ H d - £Y e L ■ H d . (2) 

After considering appropriate soft breaking terms, one can write down the Higgs potential 
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as 



- (B»H°H° + h.c.) + \{gl + gl){\Wf - \H^)\ (3) 

where m# u , TriH d and B are soft parameters with mass dimension, terms proportional to 
|/i| 2 come from the superpotential (so called F-term), and the last term comes from gauge 
symmetry (called D-term). This potential indicates that, after the electroweak symmetry 
breaking, the //-parameter is related to the Higgs vacuum expectation value (vev) and so 
its value should be around the electroweak scale. But on the other hand, since /i as the 
only parameter with mass dimension appears in the superpotential, it should naturally lie 
at the SUSY breaking scale. Such a tremendously large scale gap is usually referred as the 
u-problem Il9l . 

□ 

Different from the general MSSM where all soft breaking parameters are independent [20] , 



the CMSSM 22[ assumes following universal soft breaking parameters at SUSY breaking 



scale (usually chosen at the Great Unification scale) 

Mi/a , M , A ,tan/3 , sign(/x), (4) 

with Mi/2, Md and Aq denoting gaugino mass, scalar mass and trilinear interaction coefficient 
respectively, and evolves the four parameters down to weak scale to get all the soft breaking 
parameters of the low energy MSSM. In this sense, the parameter space of the CMSSM 
should be considered as a subset of that for the MSSM, and so is its phenomenology. 

NMSSM and nMSSM : In order to solve the /i-problem in the MSSM, various singlet 
extensions of the MSSM were proposed in history, and among them the most well known 
models include the NMSSM and nMSSM. The superpotentials of these two models are 



respectively given by [23j, l24j ] 



W NMSSM = W F + XSH U ■ H d + -S\ (5) 

3 

jynMSSM = Wf + \s H u ■ H d + £ F M*S , (6) 

where A, k and £p are dimensionless parameters of order 1, and the dimensionful parameter 
M n may be naturally fixed at weak scale in certain basic framework where the parameter is 



generated at a high loop level 24j . One attractive feature of the both models comes from the 



fact that, after the real scalar component of S develops a vev (S), an effective /x parameter 



of (9(100GeV) is naturally generated by /i e // = X(S). Another attractive feature is due to 
the Higgs self interactions, the squared mass of the SM-like Higgs boson gets an additional 
contribution X 2 v 2 sin 2 2/3, and furthermore, it can also be enhanced by the doublet- singlet 



mixing 



121 ] . This feature is very useful to alleviate the fine tuning of the models 12]. 



Noting that the tadpole term in the nMSSM only affects the masses of the Higgs bosons, 
one can easily conclude that the interactions in the nMSSM are identical to those in the 
NMSSM with k — 0. This fact enables us to modify the package NMSSMTools to study the 
phenomenology of the nMSSM 25[. As analyzed in 25|, in the nMSSM the lightest neutralino 
as the dark matter candidate is light and singlino dominated, and it must annihilate through 
exchanging a resonant light CP-odd Higgs boson to get the correct relic density. As a result, 
the SM-like Higgs boson tends to decay dominantly into light neutralinos or other light 



Higgs boson so that its total width enlarges 
different from the predictions of the NMSSM 



greatly. This character is rather peculiar and 



26]. 



Formula in calculation: In order to study the h — > Z7 decay and its correlation with 
the h — > 77 decay in SUSY, we define following normalized rates at the LHC and the 
international linear collider (ILC) as 

Zl ' " asuipp -+h^Zi) a*£ Br SM (/> Z 7 ) \C^) T^(h) T tot (h) ' 1 } 

a{pp _> h _> 77 ) atot Br(/^77) ^(c h9 X T^jh) Iggfo) 
77 " *m(pp h 77) 0% Br S M(^ -> 77) " \ Cg£J ' T™(h) ' V tat {h) ' ^ 1 

ff (e+e- -> Zfe -> ZZ 7 ) „ (ChzzY Wfe) Tf%(h) 
Zl ~ a SM (e+e- -+ Zh -+ ZZ7) ~ \CfgJ ' T s Jf(h) ' T tot (h) ' [J) 



r a(e+e -»■ gfe -» Zbb) ( C hZZ 

^bb = - 7-7— -^77 ~~ ^777 - TTsm - ' 77sI7777 ' "F 7T7- U U J 



a SM (e + e- -+ Z/i Z66) V^zz J " r ^( /i ) " r ^) 
Here C/j 55 and Cy^z are the couplings of the Higgs boson to gluons and Zs respectively, 
and Tz-y(h), T 77 (/i) and T b i(h) are the widths for the decays h — > Zj, h — > 77 and h ^ bb 
respectively. In getting these formula, we neglect SUSY radiative correction to the signals. 
These corrections are expected to be at few percent level given that heavy sparticles are 
preferred by current LHC experiments. 

In SUSY, the decays h — > Z7 and h — > 77 get new contributions from the loops medi- 
ated by charged Higgs bosons, sfermions (including stops, sbottoms and staus) as well as 
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charginos, and consequently, the formula of and T 77 are modified by 



r 77 (/i) 



G\m 2 w a m| 

64 7T 4 



128v^tt 3 



^ + ^ 7 + Af + A^ ± + A 7 ! 



(12) 



where A4 (i = W,t, f , H ± ,x ± ) denote the contribution from particle i mediated loops, and 
their explicit expressions are listed in Appendix. Note that our expressions differ from 
those presented in in two aspects. One is we have an overall minus sign for the new 
contributions Ah±, Aj and A x ±, and an additional factor 2 for the sfermion contributions. 
This sign difference was also observed recently in [21]. The other difference is we have 
included in a neat way the contributions from the loops with two particles (such as f\ and 
fi o r xt and xf) running in them. Such contributions were considered to be negligibly small 
[271 ]. but our results indicate that sometimes they may play a role. Also note that in the well- 
known SUSY package FeynHiggs [29[, the decay h — > Zj is actually not calculated, and in 
the package NMSSMTools |28J, the decay is calculated only by considering the contributions 
from the SM particles and the charged Higgs boson. We improve these packages by inserting 
our codes for h — > Z'j into them. 



III. NUMERICAL CALCULATION AND DISCUSSIONS 

In our calculation, we first scan over the parameter spaces of the four models by imposing 
various experimental constraints, then for the surviving samples we investigate the features 
of the h — > Z'j and h — > 77 decays. In performing the scans, we adopt the strategy described 
in [3], and we have following improvements 

• Noting the decay rate h — > Z~/ is sensitive to the chiral mixing angle Of in the scalar r 
sector (see below), we treat M L3 and M E3 as independent parameters and vary them 
in the ranges: 

100 GeV < M L3 , M E3 < 1 TeV. (13) 

• We use the latest results of the XENON100 experiment in dark matter direct 
detection [3^], the LHCb measurement of B s — > l^ + MZ 31 | and the CMS search for 



non-SM Higgs boson from the channel H/A — > t + t (32J to limit the parameter spaces. 
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FIG. 1: The scatter plots of the surviving samples in the four models, projected on the Rz-y — 
plane. Here i?^ 7 = a(pp — > h — >■ Z^/)/asM(pp — > h — > Zj) denotes the normalized Z"f signal rate 
in the SM-like Higgs boson production at the LHC, and the black line represents its upper limit 
set by the CMS collaboration 

The samples obtained in this way predicting 123 GeV< < 127 GeV, are able to explain 
the muon magnetic moment anomaly at 2a level and meanwhile satisfy various experimental 
constraints. Note for the NMSSM samples, we restrict A to vary from 0.5 to 0.7. The 
property of the SM-like Higgs boson for these samples, as pointed out in |12j, is quite 
different from that of the MSSM. 

In Fig.l we project the surviving samples on the plane of the signal rate at the 
LHC versus the SM-like Higgs boson mass in the four SUSY models. We also show the 



CMS bound on the signal rate in the figure 18[. This figure indicates that compared with 



its SM prediction, the Z^ signal rate in the MSSM and NMSSM may be either enhanced 
or suppressed, and in particular, the NMSSM can alter the rate by more than 50%. In 
contrast, the rate is always slightly suppressed in the CMSSM and severely suppressed in 
the nMSSM. We checked that, for the CMSSM the suppression is due to the increase of the 
partial width of h — > bb ll|, while for the nMSSM, it is due to the open up of new decay 
channels h — >• \Xi AA, where \ and A denote light neutralino and CP-odd Higgs boson 
respectively 25| . This figure also indicates that the SM-like Higgs boson mass can not reach 
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FIG. 2: Same as Fig.l, but showing the normalized sparticle contributions to the amplitude of 
h — > Z7 in the MSSM and NMSSM. The magenta bullets and blue triangles represent the samples 
with R^-y > 1 and i? 77 < 1 respectively. 

the measured central value of 125 GeV in the CMSSM. The reason is that we have required 
the theory to explain the muon magnetic moment anomaly so that the parameters M and 
M1/2 are upper bounded [11]. 

Next we focus on the MSSM and NMSSM. In Fig. 2 we exhibit the contributions of 
different sparticles to the amplitude of h — > Zj in the two models. Since the sbottoms and 
charged Higgs bosons have little effect on the amplitude, we do not show their contributions. 
This figure tells us following facts: 

(1) In the MSSM, the potentially largest contribution comes from f loops, which can alter 
the SM amplitude by about 10%. In contrast, the t contribution is small, usually 
changing the amplitude by less than 3%. This feature can be well understood by the 
formula listed in Appendix. Explicitly speaking, in order to get a significant sfermion 
contribution, one necessary condition is Y hLR sin 26 j/ m 2 ^ should be as large as possible, 
where YhLR denotes the chiral nipping coupling of higgs to sfermion, 9j is the chiral 
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FIG. 3: Same as Fig. 2, but showing the correlation of SUSY particle contribution to h — > 
channel with that to h — > 77 channel. 

mixing angle and m; denotes the lighter sfermion. As far as scalar top quark sector 
is concerned, a relatively light t\ must accompany with a heavy ti in order to predict 
rrih — 125GeV. Then even though A t in this case may be very large, the chiral mixing 
angle 9i is usually small, and consequently the t contribution can never get significantly 
enhanced. While for the f sector, both the parameters Mis and Me3 are unlimited, 
and one can choose light fs and an appropriate Of to maximize the contribution. In 
this process, the value of //tan/3 and the splitting between M L3 and M E3 play an 
important role in determining the contribution size. 

(2) In the MSSM, the chargino contribution is small and can only reach 3% and 0.5% for 
i? 77 < 1 and i? 77 > 1 cases respectively. The reason is in the MSSM, \x must be larger 
than 200GeV and 800GeV for the two cases, and consequently, the lighter chargino is 
gaugino-like and its coupling with h is small. 

(3) In the NMSSM, the potentially largest contribution to the decay comes from chargino 
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FIG. 4: Same as Fig.2, but projected on the plane of r|^ bY /r|™ 1 versus . 

loops with its size reaching 6% at most, while the magnitude of the f contribution is 
always smaller than 1%. This is because in the NMSSM with large A, \x is preferred 
to vary from lOOGeV to 250GeV and tan j3 is usually less than 10 121 ] . As a result, the 
hx^X^ coupling is enhanced by the sizable Higgsino component of x 1 * 1 , while the /it£tr 
coupling can not be pushed up by small /xtan/3. Figj2]also indicates that, due to the 
singlet component of h in the NMSSM, the hbb is suppressed and so is the partial width 
of h — > bb. Consequently, Rz-y may be quite large even when the chargino contribution 
is small. 

In Figj3l we show the correlation of the amplitudes for h — > Z7 and h — > 77. This figure 
indicates that in the two models, the t contribution to the amplitude of h — > Z7 correlates 
with that of h — > 77 by approximately linear relations, and even better in the chargino 
contribution. This is understandable since both the corrections come from the the same 
particles. This figure also indicates that the linear relation is spoiled for the f contribution 
in the MSSM. We checked that this is because the 9 T dependences of the two amplitudes 
are quite different, and meanwhile 9 T in this work varies over a broad range. For example, 
if Mlz ~ Me3, we have 9 T ~ 7r/4. In this case, the Zr*fi coupling and the /if 1 *T2 coupling 
approach zero by accidental cancelation so that the f contribution to h — > Z7 is very small, 
while the contribution to h — > 77 is maximized since it is proportional to sin2# r . On the 



10 





MSSM 

, •. *.>*. * 







- 


NMSSM 







2.0 0.0 



1.0 



2.0"" 



nMSSM 



/ CMSSM 



rr 



FIG. 5: Same as Fig.l, but showing the correlation between Rz^ = cr(pp — > h — > Z7) / 'cjsm(pp — > 
h — > Zj) and i? 77 = a(pp — > h — > 77) / o"sm(pp —> h 77) in the four models. 

other hand, if |M| 3 -M| 3 | > m T /ztan /3, which means # T — > 0, the contributions to h — > Z7 
and h — > 77 are all suppressed since their dominant contributions are proportional to sin 29 T . 

Considering i?^ 7 is mainly determined by the partial width of h — > and the total width 
of the SM-like Higgs boson, we present in FigJUthe ratio of r|^ SY /Tf 1 ^ versus the ratio of 
rttrta^/Ftotai f° r the two models. The left panel indicates that for almost all the samples 
in the MSSM the total and Z7 partial width of the SM-like Higgs boson is larger than the 
corresponding SM predictions. These features originate from the enhanced decay width of 
h — > bb and the constructive contributions of the SUSY particles to the decay h — > Z^ 
respectively. Interestingly, the largest increase of occurs when T tot ~ T tot (SM). The 
right panel indicates that, in order to enhance the Z7 signal in the NMSSM with large 
A, the SM-like Higgs boson must have a small singlet component to suppress the partial 
width of h y bb. In this case, even though T^ 7 is suppressed too, we have r^ 7 /r^ 7 (S'M) > 
Ftot/Ftot(SM). Moreover, as mentioned before, the T z -y can be significantly enhanced by 
the chargino contribution. 

Now let's summarize the Z7 and 77 signal rates at the LHC in different SUSY models. 
This is shown in Figj5j From this figure we have the following conclusions: 

(a) In the MSSM, although the partial width of h — > Z~f is enhanced by at most 20% (see 
FigHJ), due to the increase of the the Higgs width ad also due to the slight suppression 
of the hgg coupling 12], i?^ 7 > 1 is achieved only for i? 77 > 1.25. The maximal 
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FIG. 6: Same as Fig.l, but showing the correlation between /C^ 7 = a{e + e~ — > Zh — > 
ZZ7)/c7sM(e + e _ — > Zh — > ZZj) and K, b i = a(e + e~ — > Zh — > Zbb) / a$yi(e + e~ — > Zh — > Zbb) 
in the four models. 

values of Rz^ and i? 77 are 1.2 and 2 respectively. Among the sparticle contributions 
to Rz^ and _R 77 , the f loops may play the dominant role. The difference for the 
two observables come from the fact that, for 9 T ~ 7r/4, the f contribution to i? 77 is 
maximized, while that to h — >■ Z7 approaches zero. 

(b) In the NMSSM with large A, the sparticle contributions to the amplitudes of h — > Z7 
and h — > 77 are at most 10%, and the main mechanism to alter the predictions of 
i?z 7 and i? 77 is through tuning the hbb coupling by the singlet component of h. As a 
result, Rz-y and i? 77 are highly correlated, and both of them vary from 0.2 to 2. 

(c) In the CMSSM and nMSSM, i?^ 7 and i? 77 are slightly and strongly suppressed re- 
spectively. As discussed before, in the CMSSM the suppression is due to the increase 
of the partial width for h —> bb, while in the nMSSM, it is due to the open up of the 
exotic decay channels of h, i.e. h — > AA. 

Since the e + e~ collider provides a clean environment to observe the decays h — > Z7 and 
h — > bb, we also investigate their normalized rates at the ILC The corresponding results are 
shown in FigEJ This figure exhibits following features 

(a) In the MSSM, a suppressed Z7 signal (compared with its SM prediction) tend to 
correspond to an enhanced bb signal, and an enhanced Z7 signal requires the bb signal 
rates to be roughly at its SM prediction. In any case, the enhancement factor for the 
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two signals are less than 1.2. Note there exit few cases where the bb signal rate is 
slightly suppressed. 

(b) In the NMSSM with large A, the normalized bb signal rate is less than 1.1, and in some 
cases, it may be significantly suppressed. In contrast, the Z7 signal rate can be either 
greatly enhanced or severely suppressed. In the enhancement case, the bb signal rate 
is usually less than its SM prediction, and the greater enhancement corresponds to the 
stronger suppression. 

(c) In the CMSSM, both the signal rates are roughly equal to their SM predictions, and 
in the nMSSM, both the rates are strongly suppressed. 

Furthermore, we checked that the 77 signal rate at the ILC has similar dependence on the 
bb rate for the four models. 

IV. CONCLUSION 

In this work, we investigated the rare decay of the SM-like Higgs boson, h — > Zj, and 
studied its correlations with h — > 77 in low energy SUSY models. We performed a scan 
over the parameter space of each model by considering various experimental constraints and 
presented our results only in the experimentally allowed parameter space. We have following 
observations: 

(i) In SUSY models, the correction to the amplitude of h — > Z7 induced by a certain 
sparticle is correlated with that of h — > 77, and the latter is usually several times 
larger in magnitude than the former. 

(ii) In the MSSM, the net sparticle contribution to the amplitude of h — > Z7 is construc- 
tive with the corresponding SM amplitude and can enhance the corresponding SM 
prediction by at most 10%. Consequently, the Z7 signal rates at the LHC and the 
ILC can be enhanced by 20% at most. As a comparison, the 77 signal rates can be 
enhanced by a factor 2 due to large f contributions. 

(iii) In the CMSSM, due to the slightly enhanced h bb partial width, the Z7 signal rates 
at the LHC and ILC are slightly below their SM predictions, which is quite similar to 
the behaviors of the h — > 77 channel. 

13 



(iv) In the NMSSM with large A, the SUSY corrections to the amplitudes of h — > Z7 and 
h — > 77 are at most 10%, and to get significant deviation of the two rates from their 
SM values, the hbb coupling must be moderately suppressed by the singlet component 
of h. In this model, the two rates are highly correlated and vary from 0.2 to 2. 

(v) In the nMSSM, the signal rates of h — > Z7 and h 77 are both greatly suppressed 
due to the possible exotic decay h — > x°X°, A A. 

Since different SUSY models predict different correlation behaviors between the Z7 and 
77 rates, and some non-SUSY models also predict rather different correlation behaviors js], 
such a correlation may play an important role in testing the Higgs property and distinguish- 
ing the relevant new physics models in future. 
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In SUSY, the decays h — > Z7 and h — > 77 get new contributions from the loops mediated 
by charged Higgs bosons, sfermions (including stops, sbottoms and staus) and charginos, 
and as a result, the formula of Y Zl and T 77 are modified by 
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T(h Zi) 




am 



( 1 - 4) 3 <J + A Z7 + A% + Af + A z l 2 (14) 
\ m hJ " i 



2 



r(h 77) 



128V2tt 3 



A^ + Ar + A^+AV + A 1 ! . (15) 



The expressions of AJ 1 are relatively simple and are given by 



A™ = g hVV A 1 (r w ), AT = g h uN c Q 2 t A 1/2 (r t ), AJ = £ ^N c Q)A Q {r^ 



i fi 




-Ao(t h ±), A 1 ! = Y] —g + A 1/2 (T ±) 



(16) 



11 



where 73 = 4m 2 /m 2 ! , gtxY denotes the Higgs coupling with particles XY and A Q ,A 1 / 2 and 

A\ are loop functions with scalars, fermions and gauge bosons running in the loop. The 
explicit expressions of ghXY and A functions are given by 

Qhvv = Shi sin P + S h2 cos/3, (17) 

9htt = S hl / sin f3, (18) 

9hhh = 2( v / 2Gf) 1/2 ( Yfeii C ° S2 6f + YkRR Sin2 9f + YhLR 29 f ~) ' ( W) 

9 h f 2 f 2 = 2(^/2^)1/2 ( YfeLi sin2 / + YhRR cos2 / ~ Y/lLi? sin 26 '/) ' ( 20 ) 

<7wf+H- = f*>,(nj& + nf 3 ) - Utt n]g - ^njfi) + v^Xl + ^=An) 



L h3 



4V2 

ng = 2S , w C r i (7 fc , Ci = cos /?, C 2 = sin /3, 



+^[««(nH-nS) + t; i (ng-ng; 
+-4= ( «„(nj£ + n 22 + + v d (ul\ + n 22 + 2njg) 



w - = (vd^ Ql ' t;d= (72^ Cl ' W ' = " /A ' (21) 

4+xJ = ^(^i^i^2 + SnUaVp), gl xU - = ^{S hl U jX V i2 + S h2 U j2 V tl ), (22) 

A (x) = -x 2 [af 1 - fix' 1 )] , (23) 

A 1/2 (x) = 2x 2 [x- 1 + (x- 1 - l)/^" 1 )] , (24) 

A x (x) = -x 2 [2x~ 2 + 3X" 1 + 3(2z _1 - , (25) 

where S is the 2 x 2 (3 x 3) rotation matrix of MSSM (NMSSM) higgs mass matrix under 
the basis (H®, H%, S), h in Shi denotes the row index of the SM-like Higgs, Yhxy denotes 
the SM-like higgs coupling to sfermion interaction states, U, V denote the rotation matrices 
of the chargino mass matrix, and f(x) is defined by f(x) = arcsin 2 \/x. 

As for Af 1 , due to m z ^ and the existence of ZXY [X ^ Y) couplings, their expres- 
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sions are rather complex 



Aw = 9hvvCwAi(r w , X w ), Af 1 = g h ttN c Qt— A 1/2 (T t , X t ), 



,2 



A H ± = — o — 2 — WH±A)(rff±,AH±; 

zm 



' f'f 

45; =^^iv c Q / , /l2 (4 1) + 4 2) ), 

h h 

z = y ^vw _A 1/2 (t+,X+) + A z 2 . 

xf ^ m ± h ^ ^ 1 K x > x ' ' xtx 2 

xf;m,n=L,R x » 

Zl 2m w ( L t R R (i) 

y Xi x 2 v 



+ (#, L + -c/f + _ + #f + -g L 7 + _)^S 2 /l , (26) 

where Aj = 4m 2 /m| and the coupling coefficients of /i and Z are given by 

0* = 27* - 4Q t s 2 , ^ ± = (c 2 -s 2 )/c ro , (27) 

W/l = (7> 3 cos 2 B f - - Q f s 2 J/c w , v h = (Tf sin 2 9 f - - Q f s 2 J/c w , (28) 

V fl2 = (~ T f 8111 6 f COS ^/)/ C -' ( 29 ) 

= 2{ ^g f) i/2 Qofo" - Yhll) sin 29 f + cos 20;) , (30) 

4tfxr = ^ + 1 _ 2s -)/( 2c -)' sfitfxi = ^ + 1 _ 2s -)/( 2c -)' ( 31 ) 

= ^ + 1 - 2s -)/( 2c -)' 9 Zxtx - = (U 2 21 + 1 - 2 S 2 )/(2c w ), (32) 

= ( V nV2i)/(2c w ), g* xtj _ = (U ll U 21 )/(2c w ). (33) 

In above formula, we have defined some new functions as 

A®=4m h m h {C$ + Cg>) 

A« = 2m xt ^m-^[{2C^ + 3Cg } + Cf ) + (2Cff + C 

4/2 = 2 "V v/^V^F (^f + + Cf ) + (2Cg> + Cf)) (34) 

with Cy denoting three points loop functions introduced in 33[, and = 

C{P^P z ,m h ,m h ,m h ), = C^\ mf ^ mfV = ^U^^^, and C( 4 ) = 

16 



,(4). 
y 12 , 



C^\ m ,++ m , . For the special cases considered here, C« are given by 

xf xf 

Cb( P7 ,p„ mi , mi , mi ) = -/ 1 *iln|?»±f| 



(35) 



Ci 2 (p J ,p z ,m 1 ,m 1 ,m 2 ) = \ dy^—g^- (36) 



o 



C 23 {p y ,p z ,m 1 ,m 1 ,m 2 ) = J dy II — In | — - — | I (37) 



where 



b = -( m l ~ m z)( l ~ y), c= -m 2 z y(l-y) + mly + rnl(l-y). (38) 

Note that for the special case mi = m 2 , these functions can be further simplified to get their 
analytic expressions. The other functions relevant to our calculation are defined by 

A (x,y) = h(x,y), (39) 
Ai/ 2 (x,y) = h(x,y) - I 2 (x,y), (40) 
= 4(3-tan 2 ^)/ 2 (x,y)+ [(1 + 2x _1 ) tan 2 ^ - (5 + 2x -1 )] h(x,y), (41) 



with 



h{x ' y> = w^y] + w^w^ - /(!rI > | + j0w Mx ' V) - 9{r% (42) 

h{x,v) = -^L-lfix- 1 ) -/(»-')], (43) 
g(x) = V x^ 1 — 1 arcsin \fx. (44) 
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